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ABSTRACT 

We report on annual parallax and proper motion observations of H 2 O masers in 
S235AB-MIR, which is a massive young stellar object in the Perseus Arm. Using 
multi-epoch VLBI astrometry we measured a parallax of tt = 0.63 ± 0.03 mas, corre¬ 
sponding to a trigonometric distance oi D = l.SGj’jp Qg kpc, and source proper motion 
of {jLaCOsS, fis) = (0.79 ± 0.12, —2.41 ± 0.14) mas yr“^. Water masers trace a jet 
of diameter 15 au which exhibits a definite radial velocity gradient perpendicular to 
its axis. 3D maser kinematics were well modelled by a rotating cylinder with physical 
parameters: i;out = 45 ± 2 km s“^, Vrot = 22 ± 3 km s“^, z = 12 ± 2°, which are the 
outflow velocity, tangential rotation velocity and line-of-sight inclination, respectively. 
One maser feature exhibited steady acceleration which may be related to jet rotation. 
During our 15 month VLBI programme there were three ‘maser burst’ events caught 
in the act which were caused by an overlapping of masers along the line of sight. 
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1 INTRODUCTION 

The formation of massive stars has been one of the most ex¬ 
citing astronomical problems of the recent decade. Forming 
stars above 8 Mq cannot be explained by simply upscal¬ 
ing low mass star formation theories since radiati on pres¬ 
sure limits further accretion llPalla fc Stahleil[l993l l. Of the 
mainstream theories it is emerging that disk aided accretion 
is a strong candidate for explaining the formation of stars up 
to at least B type, as observational evide nce of disks associ¬ 
ated with such acc reting stars is found dHirota et al.l[20l3 : 
iBeltran et al.l[2004l l. Such protostellar disks are thought to 
be able to shield accreting material from stellar radiation 
and tra nsport it down to radii of ^100 au from the c entral 
object dZinnecker fc Yor^l2007l : ICesaroni et akllioOTh . 

However, one necessity of disc accretion theory which 
has yet to be properly confronted with observations concerns 
how angular momentum is carried away from the protostel¬ 
lar disk to enable accretion. Proposed methods of removing 
angular momentum include Alven waves in the disk mate¬ 
rial, magnetic breaking in the ionic component of the rotat¬ 
ing disk, and by physically carrying rotating material away 
from the disk vertically v ia magnetocentrifugally driven ro¬ 
tating jets (see review bn lKonigl fc Pudrit j|200(]l l. 

In addition to disk observations it is therefore also ap¬ 
pealing to look for evidence of rotating jets near massive 
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young stellar objects (MYSOs) in order to uncover their 
role in the angular momentum budget of accreting stars. A 
demonstration of su ch an investigation is the Ori-S6 MYSO 
dZapata et ahllioidl . This kind of work poses an observa¬ 
tional challenge since MYSOs are deeply embedded an d 
rarely have observable collimated jets dNavarete et al.f2015l l. 
We therefore turn to masers in our investigation of rotation 
in MYSO jets. 

Collisionally excited H 2 O masers are associated with a 
diverse range of kinematic phenomena in the young stel¬ 
lar objects (YSOs) of star forming regions (SFRs), in¬ 
cluding jets, wide-angle outflows a nd expanding shells (see 
for example iTorrelles et al.l l2012l l. The bright and com¬ 
pact masers are ideal for very long baseline interferome¬ 
try (VLBI) observations, which can penetrate the obscur¬ 
ing parent cores in whi ch massive stars form and reveal the 
motio ns of young jets dMoscadelli et al]|201ll : iFuruva et al.l 
I2OOOII . Multi-epoch observations allow determination of the 
internal 3D kinematics at shocked regions via combination of 
line-of-sight and proper motion measurement. Furthermore 
the annual parallax of a maser source provides a precise es¬ 
timation of its distance, which is crucial for evaluating its 
physical properties. As such VLBI is an indispensable tech¬ 
nique for studying the formation of stars and provides a 
mechanism for the search for rotating MYSO jets. 

VLBI observations in this paper were carried out 
using VERA (VLBI exploration of radio astrometry) 
llKobavashi et al)ll2003h which is a Japanese VLBI network 
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dedicated to maser astrometry of primarily Galactic sources. 
The bright H 2 O maser in S235AB was chosen as a VERA 
targ et for its interest ing ‘maser burst’ behaviour (up to 120 
Jy) dFelli et al.ll2007l l which has a dominant blue-shift maser 
(DBSM) spectrum. Its annual parallax was also sought for 
it’s contribution to the mapping of the Perseus Arm. 

In addition to H 2 O masers, th e S235AB region 
also exhibits class I methanol masers dKurtz et al.l l2004l : 
IPratap et al.l l2008l l and two bipolar outflows of HCO^, 
drive n from a dusty molecular c o re see n in 450 and 850 
fim (iFelli et al.l l2004h . iFelli et ^ (l200di proposed a class 
I mid infrared source, S235AB-MIR, as the outflow pro¬ 
genitor from Spitzer images. By SEP mode lli ng of sub¬ 
millimeter continuum data from Belli _gt__al, ([2004 ) and 
Spitze r mid-infrared photometries. fPewangan fc Anandaraol 
(l201lll showed S235AB-MIR to be the only massive YSO in 
the S235AB cluster, with a stellar mass of 11 Mq and still 
accreting. 

N o centimeter emiss i on has been dete cted in S235AB- 
MIR (iTofani et al.l 1 19951 : IFelli et all I2OO6II , suggesting no 
significant developmen t of a n ultra compact Hu region 
fUCHIIRl. ISaito et al.l (I2OO7II estimated the most massive 
member of the S235AB cluster to be of spectral class B1 
from C^®0 observations. Thus, a plausible interpretation 
above picture is that of a deeply embedded accreting B-type 
star which is too young to produce a detectable UCHIIR. 


2 OBSERVATIONS AND DATA REDUCTION 

Simultaneous VERA observations of S235AB-MIR and the 
continuum re ference source J0533-b3 451 were made in dual¬ 
beam mode l|Kawaguchi et ai]l200d l. Observing the maser 
and reference source simultaneously removes the need to 
slew antennae and interpolate phase measurements between 
sources, enabling excellent determi nation of the dyn amic 
troposphere phase contribution fsee lAsaki et al.ll2007l l. 

Phase tracking centers were set to (a, (I)j 2000.0 = 
(05‘’40“53".38445496, -b35°41’48”.4470) and (q,( 5)j2ooo.o = 
(05‘’33™12=.7651060, -b34°51’30”.336990) for S235AB-MIR 
and J0533-I-3451, respectively. The continuum source 
J0533-I-3451 is listed in the VLBA calibrator list dPetrovI 
l2012f) and was observed at an unresolved A-band flux of 
~20 mjy in our observations with VERA. Intermittent ob¬ 
servations of BL Lac, DA55 or 3C84 were made every 1.5 hrs 
for bandpass and group delay calibration. Typical observing 
sessions were 8 hrs long, providing 2.5 hrs on-source time 
and sufficient uti-coverage. 

Left-handed circular polarisation signals were sampled 
at 2-bit q uantisation, and fi ltered with the VERA digital 
filter unit dlguchi et al.ll2005l L The total bandwidth of 256 
MHz was divided into 16 intermediate frequency (IF) chan¬ 
nels, each with a bandwidth of 16 MHz. One IF was allocated 
to the maser signal, assuming a rest frequency of 22.235080 
GHz. The other 15 IPs, in adjoining frequency, were allo¬ 
cated to J0533-I-3451. Signal correlation was carri ed out us¬ 
ing the Mitaka EX correlator (IChikada et ahlfloOll L Channel 
spacings of 31.25 kHz were used for the maser data, corre¬ 
sponding to a velocity resolution of 0.42 km s“^. 

Global positioning system (GPS) measurements of the 
atmospheric water vapour content at each station were 
used to refine an a-priori model of expected atmospheric 


Table 1. Summary of observations made with VERA. 


Observation 

Epoch 

Date 

MJD 

Detected 

spots 

1 

2013 Jan 29 

56321 

17 

2 

2013 Feb 21 

56344 

18 

3 

2013 Mar 15 

56366 

26 

4 

2013 Apr 21 

56403 

12 

5 

2013 Sep 29 

56564 

30 

6 

2013 Nov 02 

56598 

24 

7 

2013 Dec 02 

56628 

26 

8 

2014 Jan 28 

56685 

21 

tg 

2014 Mar 11 

56727 

19 

tfO 

2014 Apr 29 

56776 

19 


f These epochs were not used in parallax determination. 


delay llHonma et al.l l2008l f. Solutions were applied post¬ 
correlation. 

VERA observations were conducted in 10 epochs, 
closely spaced to support maser identification in view of the 
extreme variability of emission. The observation calendar is 
summarised in Table [T] An error (possibly caused by a shift 
in the adopted antenna positions) in our correlation settings 
in epochs 9 and 10 resulted in inaccurate astrometry such 
that data from these epochs could not be used in measure¬ 
ment of the annual parallax. These data were only used in 
maser distribution and line of sight velocity analysis. 

All data were reduced using the Astronomical Image 
Processing System (AIPS) developed by the National Radio 
Astronomy Observatory (NRAO). 

Data were reduced using the inverse phase- referencing 
metho d for VERA data, which was introduced in llmai et al.l 
1 I 2 OI 2 II . In this approach the bright maser emission is used 
to calibrate dynamic tropospheric fluctuations in phase and 
rate - thus acting as the phase reference - while group de¬ 
lay is solved using a bright continuum source. These solu¬ 
tions are applied to the data of a nearby continuum source 
whose precise co-ordinates are already known. The relative 
position of the maser to the nearby continuum source gives 
the accurate position of the maser. A detailed explanation, 
a flowchart and an internet link to a data reduction walk¬ 
through are given in the Appendix of this paper. 

The main advantage of the inverse phase-referencing 
technique is its ability to detect weak continuum sources 
(~10 mJv. llmai et al.ll2012fl . Furthermore, since the phase- 
rate solutions applied to the nearby reference source are also 
used to self-calibrate the maser data, the phase-referenced 
maser maps have high dynamic range. 

Emission peaks in the maser maps were ima ged using 
the CLEAN procedure, based on iHoeboml (Il974l') and de¬ 
tections were registered at a signal-to-noise cutoff of 5. Fol¬ 
lowing common nomenclature, a maser ‘spot’ refers to an 
individual maser brightness peak, imaged in a specific ve¬ 
locity channel, and a maser ‘feature’ refers to a group of 
spots which are considered to emanate from the same phys¬ 
ical maser cloud. Maser spots are categorised into features 
when they are grouped within 1 mas of another spot, and are 
continuous in velocity. Features and their associated spots 
are denoted by letters and numbers respectively in column 
1 of Table H 
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3 RESULTS 

3.1 Maser detections 

In total, 19 maser features were catalogued comprising some 
~60 spots, all of which showed notable time variability. 
Maser detections are summarised in Tabled Features can be 
separated into three main groups by comparing local stan¬ 
dard of rest (LSR) velocities to that of the molecular core, 
which has a velocity of —17 km s“^ dFelli et al.ll2004 l. The 
groups are: ‘red lobe’, ‘core velocity’ and ‘blue lobe’ masers, 
which were observed from -|-2 to -1-14, around —18, and from 
—55 to —65 km s“^, respectively. 

Red lobe and core velocity masers were detected spo¬ 
radically, and were weak - with typical fluxes of ^ 1 Jy. The 
red lobe group had only two individual maser features and 
the core velocity masers had only one maser feature. Masers 
from these groups reside loosely together in a region of about 
120 X 120 au. 

Blue lobe maser emission was strong, extremely variable 
and exhibited a ‘burst’ nature where the cross-power flux 
rose from 10° Jy to 10^ Jy on three occasions. These masers 
are blueshifted > 40 km s“^ with respect to the molecular 
core. Blue lobe masers reside in a compact region of about 15 
X 15 au at a separation of 2” from the other maser groups. 

Relative intensities and time variation of the aforemen¬ 
tioned velocity components are shown in Figure[T] The spec¬ 
trum of S235AB-MIR exhibits a classical DBSM profile, pro¬ 
viding e arly evidence that the water masers may be drive n 
by a jet (ICaswell &: Philli^l2008l : [Motogi et al.ll2013l . l2015ll . 
Figure[2]shows the velocities and distributions of all detected 
maser features, revealing a jet. A dashed line drawn through 
the central positions of the main maser groups traces a po¬ 
sition angle of PA = 217° which delineates the geometric 
axis of the system. The right-hand panel of Figure [2] shows 
a position-velocity (p-v) diagram across this axis. 


3.2 Parallax and proper motions 

Astrometric motions of masers are a combination of a linear 
component corresponding to the proper motion of the maser 
on the sky, and a sinusoidal component due to annual par¬ 
allax. These two components are separated by simultaneous 
fitting of a linear and sinusoidal function to the observed 
data. Long temporal baselines and dense sampling are vital 
to correctly separate the two components, therefore we only 
fit parallaxes using masers associated with the persistent 
features; A and B. 

We performed data fitting using 3 different approaches: 
individual, group and feature fitting. ‘Individual fitting’ 
refers to fitting the parallax and proper motion of each indi¬ 
vidual maser spot independently. Fitting is done using the 
same velocity channel at each epoch i.e. we assume non¬ 
acceleration in the line of sight. We use only the most stable 
masers that were identifiable in at least 5 epochs and persist 
for at least 1 year. 

‘Group fitting’ involves fitting the astrometry of multi¬ 
ple spots together by assuming a common distance. In this 
procedure we use all masers that are identifiable in at least 4 
epochs and are associated with features A or B. We perform 
‘group fitting’ on masers in features A and B together, and 
also separately for comparison. 



Figure 1. Scalar averaged spectra of the H 2 O emission for each 
epoch, limited to 2 Jy to show weaker detections. Lines of constant 
velocity indicate the red lobe, the core velocity and the blue lobe 
maser velocities in descending velocity order. 


‘Feature fitting’ is an approach where we determine the 
nominal position of the entire maser feature by calculat¬ 
ing the flux-weighted average position of the three brightest 
maser spots. 

The results of these three fitting procedures are sum¬ 
marised in Table [ 2 ] where astrometric positions in columns 
4 and 5 are quoted as offsets from the phase tracking cen¬ 
ter of the maser beam which was set to (a, 5)j2000.0 = 
(05‘’40“53".38445496, -b35°41’48” .4470). 

‘Individual fitting’ was done on 3 masers in feature A. 
Each maser spot gave consistent parallax and proper motion 
values. The ‘group fitting’ of maser spots in feature A gave 
results that were consistent with the average of the ‘indi¬ 
vidual fitting’ results. However, when using ‘feature fitting’ 
the same consistency could not be obtained. For thermal 
sources the flux of each component should be steady, thus 
the nominal position of the feature should trace the gravita¬ 
tional center. However, feature A exhibits multiple velocity 
components as can be inferred from column 3 of Table [2] (see 
table footnote) and any changes in the relative brightnesses 
of components will affect the determination of the nominal 
astrometric position of the feature. Therefore, regarding fea¬ 
ture A, fitting is best done using spots rather than features. 

Multiple maser components also exist in feature B as 
is confirmed by its maser burst nature (see subsection 4-^)- 
As such, ‘feature fitting’ was also inappropriate for feature 
B. Emission was not detected for a whole year in any one 
channel barring ‘individual fitting’, however ‘group fitting’ 
for this feature gave a parallax that matched very well to 
the results for feature A. 

Maser emission in feature C per individual epoch is nar¬ 
rower in velocity than that of features A and B, and does 
not exhibit multiple velocity peaks or a burst nature. This 
suggests that feature C has a simpler physical and velocity 
structure than features A and B. However, feature C is seen 
to increase in velocity over time (see subsection 4-5)- As 
a result, the temporal baseline of emission seen in any one 
maser channel was too short for astrometry using ‘individual 
fitting’ or ‘group fitting’. The proper motion of this maser 
feature was obtained through ‘feature fitting’, its parallax 
value was consistent with the aforementioned results from 
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features A and B, however the parallax was not used as it 
violates the criterion that masers must not be accelerating. 

In order to conclude the annual parallax of S235AB- 
MIR we decided to use ‘group fitting’ fo r maser spots in 
featu res A and B together. Error floors fsee lHachisuka et al.l 
I 2 OO 9 I I of 0.19 and 0.15 mas were added to RA. and Dec. 
positions, respectively, in order to achieve a ~ ^ fit- The 
annual parallax of S235AB-MIR was tt = 0.639 ± 0.033 mas 
corresponding to a trigonometric distance of D = 1.564 q q 8 
kpc, firmly placing it in the Perseus Arm. 

The failure of ‘feature fitting’ for features A and B 
demonstrates that the type of fitting procedure used can 
affect parallax estimates, especially in the presence of unre¬ 
solved structure. Proper motions shown in Table [5] are those 
determined in by group fitting, assuming this common dis¬ 
tance. Maser astrometries and their components are shown 
in Figured 


3.3 Internal and systemic motions 

The observed proper motions of Galactic masers are a com¬ 
bination of internal motions inherent in the SFR, and the 
systemic motion of the SFR with respect to the Sun. 

Considering Galactic motion, S235AB-MIR resides 
in the Galactic anticenter direction {I = 173.7°), where 
systemic proper motions ar e smal l. Using the Galactic rota¬ 
tion curve from iReid et al.l ll20l4 l: 0 = 0o — 0.2(i? — Ro), 
where Ro = 8.34 ± 0.16 kpc, 0o = 240 ± 8 km s“^, and 
using our distance to S235AB-MIR; 1.56 kpc, we calculate 
0 = 240 km s“^. From this, we estimate the proper motion 
cont ribution due to Ga lactic rotation using a formulation 
from lAndo et ahl (1201 il l: 

fj,icosb — ~ i^')^ocos I — 

Which gives ^iCos6 = —0.019 mas yr~^. We must also cor¬ 
rect for the motion of the Sun with respect to the LSR, which 
is taken as ( U, V, W )(o = (-1-10. 3,-1-15.3,-b7.7) km s~^ from 
lAndo et al.l (l201ll ~) (based on iKerr fc Lvnden-Belll Il986ll . 
These contribute an apparent motion of (uaCosS, fis) = 
(-1-0.27,—2.41) mas yr~^. By subtracting this motion from 
the observed proper motions we isolated the internal mo¬ 
tions of water masers in S235AB-MIR which are shown in 
Figure [21 

Alternatively we sought to infer the internal proper mo¬ 
tions by use of a reference maser. Feature F has a velocity 
that is well removed from those of the blue and red lobe 
masers and instead has a velocity close to that of the molec¬ 
ular core; this feature seems unrelated to the jet. By desig¬ 
nating feature F as the reference maser we can separate the 
internal motions (traced by blue and red lobe masers) from 
the source systemic motion (traced by feature F). 

Feature F has a proper motion of {fiaCosS, fig) = 
(-1-0.79 ± 0.12, —2.41 ±0.14) mas yr“^ which is almost iden¬ 
tical to the correction factor calculated above from the 
LSR-frame motion. As such both the reference maser ap¬ 
proach and LSR-frame approach produce very similar in¬ 
ternal motions. In the latter approach the source systemic 
proper motion would be equivalent to that of feature F; 
{^aCosS, us) = (±0.79 ± 0.12, —2.41 ± 0.14) mas yr“^. This 


value is consistent with o ther Perseus Arm sources cata¬ 
logued in iReid et al.l (l2014l L 

Resultant internal proper motions are typically smaller 
than the line of sight velocities and motions do not ap¬ 
pear well collimated. This can be explained by projec¬ 
tion effects from a pole-on geometry, as was seen in the 
case of G353.273±0.641 which is also a DBSM jet source 
dMotogi et ahllioish . 


4 DISCUSSION 

4.1 Driving source: S235AB-MIR 

iFelli et H] (I 2 OO 6 II proposed the class I MYSO S235AB-MIR 
as the progenitor of the molecular outflows and maser emis¬ 
sion in this region. Our observations confirm that the blue- 
shifted and red-shifted velocity components, which are sep¬ 
arated by 2”, lie either side of S235AB-MIR and in agree¬ 
ment with the corresponding blue and red-shift ed lobes of 
the NNW-SSE molecular outflows seen in HCO’^ dFelli et al.1 
l2004h . Less extended blue and red-s hifted lobes , also orien- 
tated NNW-SSE, are seen in C^'^S. iFelli et al.1 d2006lj con¬ 
sider two possible interpretations of the C^'^S velocity lobes; 
a molecular outflow and a rotating disk. Positions, line of 
sight velocities, and proper motions of water masers pre- 
sented in t his w ork support the outflow interpretation of 
IFelli et H] d2006l L such that the G®^^ S gas is a component of 
the NNW-SSE outflow. Furthermore, the proper motions of 
water masers oppose the radio and millimeter source VLA- 
1/Ml as progenitor of the NNW-SSE outflow, rather sup¬ 
porting the view that the outflows are driven by S235AB- 
MIR. 

The HGO’'" emission has a distinct quadrupole morphol¬ 
ogy - one bipolar set of lobes orientates to NNW-SSE with 
the other set almost perpendicular. At the current best res¬ 
olution the outflows appear concentric. Regarding establish¬ 
ment, one possible explanation is that of a sudden change 
in the outflow direction of a single bipolar outflow driven 
by the MYSO in S235A B-MIR, due to an encounter with 
another cluster member dBallv fc Zinneck^l2005l L Another 
possibility is that of multiple bipolar outflow progenitors 
in the S235AB-MIR core, at angular sizes finer than those 
probed in the Spitzer images. Considering this scenario, we 
cannot rule out the possibility that the high velocity wa¬ 
ter masers in S235AB-MIR are driven by a companion low 
mass YSO. However, the emission maps of IFelli et al.1 d2004h 
show that the orientation and velocity structure of the G®'*S 
(dense gas tracer) emission correlates with that of the NNW- 
SSE HCO^ outflow and water masers (this work), whereas 
no G^^^S gas was seen to trace the NE-SW outflow. So, it 
is reasonable to assume that the former may therefore be 
energetically dominant and would thus likely be attributed 
to the MYSO. 
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Figure 2. Positions and internal motions (corrected for systemic motion) of masers in S235AB-MIR. Coordinates are given respective to 
the phase tracking center of the maser data. Middle shows the full scale maser distributions including the three maser groups discussed 
in the main text. A short dashed line traces the jet axis. The proper motion of feature F {vlsr = —17.91 km s“^) is too small to see. 
Left shows a zoom of the cluster of blue lobe masers. Right shows a p-v diagram of the blue lobe masers with respect to the jet axis, 
where vertical bars are the half-widths of typical maser velocity features, larger bars at offsets of H-l and -2 mas correspond to features 
A and B, respectively, which have wider velocity profiles. The thick grey line is the least squares fit and traces a gradient of 1.2 km s~^ 
mas”^. 



Figure 3. Parallax and proper motion fitting of spots (coloured arbitrarily) associated with maser features A and B in S235AB-MIR. 
Left: Sky-plane motions of the maser spots with respect to the phase tracking center. Middle: arbitrary scale astrometric offsets in Dec. 
(above and middle for Features B and A, respectively) and R.A. (below) as a function of time. Right: parallactic motion in R.A. (above) 
and Dec. (below) for all spots in features A and B after subtraction of linear proper motions. Small horizontal offsets (within ±7 days) 
were introduced to data points in the middle and right-hand plots for readability. 
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Table 2. The general properties of H2O masers in S235AB-MIR detected with VERA. 


Maser 

l^LSR 

Detected 

Ad cos 5 

A6 

flex. COS 5 


TT 

ID 

(km ) 

epochs 

(mas) 

(mas) 

(mas yr“^) 

(mas yr~ ^) 

(mas) 

A1 

-60.06 

’^*345’^*8*’^ 

-6.95 

32.45 

0.62 ± 0.27 

-1.95 ± 0.22 


A2 

-60.48 

*^3*5**g^* 

-6.92 

32.47 




A3 

-60.90 

1*3456*** 10 

-6.76 

32.88 

1.30 ± 0.29 

-2.36 ± 0.23 


A4 

-61.32 

123456789 10 

-6.72 

32.88 

1.17 ± 0.18 

-2.42 ± 0.15 

0.713 ± 0.092 

A5 

-61.74 

123*5678* 10 

-6.68 

32.91 

1.34 ± 0.19 

-2.48 ± 0.15 

0.675 ± 0.061 

A6 

-62.16 

1***56789 10 

-6.69 

33.01 

1.20 ± 0.25 

-2.43 ± 0.20 

0.558 ± 0.083 

A7 

-62.58 

12**56*89 10 

-6.65 

33.09 

1.24 ± 0.22 

-2.43 ± 0.18 


AS 

-63.00 

12**5**89 10 

-6.61 

33.09 

1.08 ± 0.23 

-2.38 ± 0.18 


A9 

-63.42 

12**5**89 10 

-6.54 

33.02 

0.95 ± 0.23 

-2.37 ± 0.18 


AlO 

-63.84 

123*5***9 10 

-6.44 

32.97 

0.60 ± 0.37 

-2.36 ± 0.29 


All 

-64.26 

*23******* 

-6.57 

32.81 




A12 

-64.68 


-6.49 

32.70 




A13 

-65.10 


-6.44 

32.74 




A14 

-65.52 

**3******* 

-6.43 

32.76 










Average 

TT = 0.649 ± 0.144 


Feature A group fitting 




1.06 ± 0.08 

-2.35 ± 0.07 

TT - 0.656 ± 0.038 


Feature fitting result 




0.36 ± 0.19 

-2.45 ± 0.10 

TT - 1.018 ± 0.093 

B1 

-60.06 


-6.33 

36.85 




B2 

-60.48 

***4*6**** 

-6.27 

36.81 




B3 

-60.90 

**34*67*** 

-6.67 

37.18 

1.26 ± 0.31 

1.47 ± 0.24 


BA 

-61.32 

*234567*** 

-6.56 

37.05 

1.34 ± 0.26 

1.10 ± 0.20 


B5 

-61.74 

*234567*** 

-6.53 

37.05 

1.54 ± 0.26 

1.17 ± 0.20 


B6 

-62.16 

123*567*** 

-6.42 

37.28 

1.43 ± 0.23 

0.86 ± 0.18 


B7 

-62.58 

123*567*** 

-6.47 

37.14 

1.45 ± 0.23 

0.92 ± 0.18 


B8 

-63.00 

*23*5678** 

-6.24 

36.98 

1.52 ± 0.22 

0.68 ± 0.18 


B9 

-63.42 

*23*5678** 

-6.15 

36.86 

1.40 ± 0.22 

0.83 ± 0.18 


BIO 

-63.84 

**3*5678** 

-6.14 

37.02 

1.42 ± 0.29 

0.63 ± 0.23 


Bll 

-64.26 

**3*5678** 

-6.10 

36.99 

1.43 ± 0.29 

0.64 ± 0.23 


B12 

-64.68 

**3*567*** 

-6.11 

36.98 

1.31 ± 0.35 

0.88 ± 0.28 


B13 

-65.10 

**3*567*** 

-6.21 

37.06 

2.17 ± 0.35 

0.64 ± 0.28 


B14 

-65.52 

****567*** 

-3.72 

37.33 




B15 

-65.94 

****567*** 

-3.99 

37.20 




B16 

-66.36 

****5*7*** 

-4.06 

37.16 





Fearure B group fitting 




1.48 ± 0.08 

0.89 ± 0.07 

TT - 0.628 ± 0.057 


Feature fitting result 




1.89 ± 0.14 

0.74 ± 0.21 

77 = 0.449 ± 0.069 


Collective group fitting of 








spots in features A and B 






77 - 0.639 ± 0.033 

Cl 

-59.22 

123******* 

-8.103 

33.71 




C2 

-59.64 

1 ********* 

-8.063 

33.65 




C3 

-60.06 

1*34****9* 

-7.935 

33.594 




C4 

-60.48 

***4****g* 

-7.568 

33.221 




C5 

-60.90 


-6.271 

33.644 




C6 

-61.32 

****56789 10 

-6.206 

33.644 




C7 

-61.74 

****56789 10 

-6.291 

33.709 




C8 

-62.16 

****56789 10 

-6.316 

33.714 




C9 

-62.58 

****56789 10 

-6.061 

33.509 




CIO 

-63.00 

******789 10 

-6.498 

33.529 




Cll 

-63.42 

******789 10 

-6.493 

33.514 




C12 

-63.84 

******789 10 

-6.528 

33.509 




C13 

-64.68 


-6.393 

33.251 





Feature fitting result 




1.37 ± 0.19 

-0.44 ± 0.15 

77 - 0.632 ± 0.077 


D 

+ 14.12 

12******** 

-1213 

1563 




E 

+2.73 

*********10 

-1155 

1599 




F 

-17.91 

*******g9 10 

-1136 

1619 

0.79 ± 0.12 

-2.41 ± 0.14 


C 

-61.32 

1234****** 

-8.19 

36.35 

1.47 ± 0.20 

-1.91 ± 0.20 


H 

-60.91 

****5*7g** 

-7.19 

36.29 

0.84 ± 0.22 

1.44 ± 0.30 


I 

-65.94 

****567*** 

-1.68 

36.30 

3.46 ± 0.20 

-2.02 ± 0.21 


J 

-55.43 

123******* 

-9.36 

32.22 

2.60 ± 0.18 

-6.18 ± 0.20 


K 

-56.28 

123*5***** 

-5.18 

26.21 

3.77 ± 0.18 

-5.61 ± 0.20 


L 

-56.70 

***4****** 

-4.38 

25.25 




M 

-57.54 


-11.14 

37.03 




N 

-61.32 

****5***** 

-1.41 

30.14 




0 

-63.86 

*********10 

-3.12 

29.39 




P 

-64.25 

****5***** 

-5.69 

39.10 




Q 

-64.70 


-5.06 

33.18 




R 

-64.70 

*********10 

-1.67 

35.90 




S 

-64.70 

*********10 

-2.16 

33.21 





Column (3): Epoch numbers indicate detection, while asterisk represents non-detection. For each epoch the brightest velocity compo- 
nent(s) in features A, B and C are indicated by underlined boldface digits. 
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Figure 4. Scalar averaged spectra of the blue lobe H 2 O emission 
at full scale, highlighting the variable burst nature. 


4.2 Maser burst 


A sudden increase i n mase r brightness is known as a ‘burst’. 
IPeguchi fc WatsonI (llQSflI 'l proposed a maser burst model 
in which an increase in path length, achieved by the over¬ 
lapping of two masing clouds along the line of sight, can 
produce an increase of two or more orders of magnitude 
in observed brigh t ness. Multi-epoch VLBI observations by 
IShimoikura et ahl (l200Rh found observational evidence of 
this process in the case of the water masers of Orion KL, 
however there are few other examples of a maser burst 
caught ‘in the act’. There were 2 strong maser bursts during 
our VLBI monitoring observations of S235AB-MIR in which 
fluxes of ~100 Jy were reached, as shown in Figure ID 

Column 3 of Table[5]shows the detected velocity compo¬ 
nents (rows) for each epoch (columns), where the brightest 
peaks in in each epoch are shown as underlined boldface 
digits for spots associated with features A, B and C. From 
the peaks we can see that feature B has a complex veloc¬ 
ity structure; velocity drifts and multiple peaks. Two peaks 
can be seen in epoch 5 (MJD 56564) at —64.68 km s“^ and 
—62.58 km s“^ which respectively correspond to the blue 
and green components seen in Figure O where feature B 
is at map offset (0,0). In epoch 6 (34 days later) the to¬ 
tal flux of feature B raised from 7 Jy to 80 Jy as the two 
components overlapped both in space and velocity. There¬ 
fore we observe t hat the bursting mechanism in the masers 
of S235AB- MIR dFelli et al1l2007 l is the same as that which 
was seen in IShimoikura et al. (1200^ . Figure [5] illustrates the 
‘pre-burst’ stage of this mechanism. 

Bursting activity was also seen in feature B at epoch 3 
in which it reached a flux of 120 Jy. The bursting mechanism 
in epoch 3 was the same as described above. 

Feature A also has a complex velocity structure. In this 
case feature A appears as a single peaked maser feature in 
epoch 1, which gradually separates in velocity into two de¬ 
fined peaks by epoch 3, and then recombines at epoch 5 - 
where a modest maser burst raised its flux to 20 Jy from its 
usual flux of 2-^3 Jy. 



Figure 5. First moment (velocity) map of water maser emission 
in epoch 5 (MJD 56564) where positions are relative to the ref¬ 
erence maser in feature B. The three main features are labelled. 
Two velocity components in feature B, at map position (0,0), are 
the progenitors to the maser burst seen in epoch 6 (reaching 80 
Jy). Contours are 3, 5, 10, 20, 40, 80 times the image rms noise 
of 200 mjy. 


The peak velocities of the overlapping maser compo¬ 
nents involved in the burst often differed by a few km s“^, i.e. 
the burst occurred at the intermittent velocity between the 
peaks. This suggests the existence of a significant amount 
of population inverted gas at velocities away from the origi¬ 
nal peaks. Such behaviour indicates that the physical maser 
cloud may have a spatially unresolved structure and pos¬ 
sibly some internal motions since the bursting components 
were part of the same maser feature. Such scales cannot 
be investigated directly at the current angular resolution of 
terrestrial VLBI arrays, however the presence of such sub¬ 
structures has recently been observed in the results of space 
VLBI observations (RadioAstron newsletter; #19 - released 
6th March 2013, and #27 - released 31st March 2015). 


4.3 H2O masers: the jet and the velocity gradient 

H 2 O masers can be collisionally excited on the surface of 
shock fronts created when a prot ostellar jet pushes into am- 
bient gas {ex: coni cal jet surfoce iMosc adelli et^ I I2000I. U- 
shaped ‘micro- 7 'ef ’ iFuruva et ahl I2OOOI) . The large velocity 
difference in blue (—66 km s“^) and red lobe (-1-14 km s“^) 
masers in S235AB-MIR indicate that these masers are likely 
tracing a jet. Elongated red masers parallel to the jet axis, 
shown in Figure [5] may represent shocked gas on the inner 
surface of the jet. 

The jet axis determined by H 2 O masers in Figure[2]is in 
good agreement with the orientation s of large s c ale H CO'*' 
and smaller scale C^'*S outflows from iFelli et al.1 ll2004h . We 
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made a p-v diagram using offsets from the jet axis plotted 
against line-of-sight velocity in the right-hand panel of Fig¬ 
ure [5] Maser features are well fit to a straight line as would 
be expected from solid body rotation. Its gradient of 1.2 km 
s“^ mas“^ corresponds to 0.77 km s“^ au~^ at 1.56 kpc. 

Velocity gradients across jet axes have been re- 
ported in several YSO observations and simulations 
JZaDateet aL 1120101. 1 


I: ICoffev et al.l [20^: IWoitas et al.l 


I2OO5I : Bacciotti et al.l ' 20021 : ICeraueira et ajT[2006ll . To ex¬ 

plain these velocity structures there is current debate 
between supporter s of t he jet rotation i nterpretation 


between supporter s ot t ne .let rotation i nterpretation 
lZapata^_et_ajj^ 2010jj_ 20151: jCoffev et ^ l2004l : IWoitas et al.l 
I2OO5I : iBacciotti et al.l 12 0021) a nd those who propose alter- 
native explanations llSokerll2005l : ICeraneira et al.ll2006lf . Al¬ 
though we believe S235AB-MIR possibly evidences a rotat¬ 
ing jet we consider other commonly suggested explanations 
of velocity gradients in jets, before discussing a model of jet 
rotation. 


4-3.1 Interpreting the velocity gradient 

Precessing jets are occasi onally found in MYSOs 
dCaratti o Garatti et al.l I 2 OI 5 I I. In this scenario the jet 
is comprised of a continuous ballistic ejection whose di¬ 
rection is ti me-dependent, causing the well-known ‘wiggle’ 
morphology. ICeraneira et al.l (l2006l l show that the velocity 
field created by a precessing jet and a rotating jet are 
difficult to distinguish in radial velocity profiles of emission 
line observations. Such conclusions are based on large 
scale gas motions encompassed within the precession angle; 
scales larger than the radius of the jet. The H 2 O masers in 
S235AB-MIR trace scales comparable to the physical radiu s 
of the jet (~20 an, iFuruva et al.|[200^ : ICoffev et al.ll2003l . 
A velocity gradient on a jet radius scale can only be formed 
if there is motion within the jet itself - thus rejecting a 
precession scenario. 

Velocity structures may also arise from dispersion at the 
tip of a jet driven bow shock. In this scenario the shocked 
gas at the tip of the je t is forced edgew ays as it decelerates 
into ambient material. iLee et akl (120001 1 modelled jet-driven 
bow shocks and found their p-v diagrams to be symmetric 
across the central axis (their Figure 26). This contrasts with 
the p-v diagram of water masers in S235AB-MIR which is 
distinctly asymmetric (Figure [2}. Furthe rmore, rad i al ve - 
locity gradients seen in the bow shocks of iLee et al. I (I2OOOII 
are several orders of magnitude lower than that observed in 
S235AB-MIR (roughly 2 x 10“® km s“^ au“^ for the case of 
HH 111 in their Figure 27). The bow shock interpretation 
is therefore inconsistent with the velocity gradient seen in 
S235AB-MIR 

If masers were produced in the shock front of an expand¬ 
ing arc it would be possible to reproduce a velocity gradient 
across the direction of collective motion if the arc were flat¬ 
tened on a plane, inclined to the observer and fanning out 
at constant velocity from a common driving source. How¬ 
ever a fanned out jet scenario is a strictly ballistic model, 
requiring the proper motions of maser features to interpo¬ 
late to a common kinematic center. In contrast, blue lobe 
masers in S235AB-MIR form a cluster, as opposed to an arc, 
and proper motions do not interpolate; some of the maser 
features in S235AB-MIR have proper motions that are par- 
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Figure 6. Plots of X^ed fitting parameters at inclination slices 
of i = 2,12, 22° using model fitting condition 1. 



allel to the jet axis while some are near perpendicular to it 
- which is difficult to explain by purely ballistic motion. 

We also note that b lue lobe masers have been known 
since lHenkel et ^ lll986ll yet the maser group remains com¬ 
pact despite differences in line of sight velocities and di¬ 
vergent proper motions. This is difficult to reconcile using 
ballistic arguments since masers would be expected to have 
dispersed more by the present day. 

Since only a modest number of maser features were de¬ 
tected in S235AB-MIR we cannot fully dismiss all ballistic 
interpretations of the observed velocity gradient using the 
current data alone and require confirmation from observa¬ 
tions that sample jet gas directly. Nevertheless, without a 
satisfactory ballistic explanation for the kinematics of water 
masers in S235AB-MIR we proceed with a rotating jet as 
our working hypotheses. As we will see, a rotating jet inter¬ 
pretation forms consistencies with the aforementioned ar¬ 
guments that had conflicted with ballistic explanations. To 
our present knowledge S235AB-MIR is the first candidate 
rotating jet traced by water masers and, if confirmed with 
further observations, presents an ideal laboratory to study 
magnetocentrifugal jets in a class I YSO with an accurate 
distance. 

We created a kinematic model of a rotating cylinder to 
compare with the sky-plane distributions and 3D (sky-plane 
and line of sight) velocities of the water masers in S235AB- 
MIR. A cylindrical surface model was chosen for its analogy 
to an interface between a jet and the ambient gas - where 
masers are excited. Also the choice of a cylinder over a con¬ 
ical surface circumvents the ambiguity of the driving source 
position (which should not matter since jets are able to de- 
collimate and re-collimate). The fixed physical parameters 
of the model were the jet position angle, PA, and radius, r, 
which were readily determined from the maser maps (Fig¬ 
ure 0 as PA = 217°, r = 7.5 au. The free parameters were 
the linear jet velocity, Uout, the tangential rotation velocity, 
Urot, and the inclination to the line of sight, i. For compari¬ 
son with observational data the modelled 3D motions were 
projected into the sky-plane and line of sight. 

Modelled values were compared with the maser data 
with regards to: positional offset from the jet axis, line- 
of-sight velocity with respect to the parent cloud ( — 17 km 
s“^) and, when present, R.A and Dec. internal proper mo¬ 
tions. We employed a ‘reduced x^’ fitting approach in which 
Xred = where N is the number of degrees of freedom, 

evaluated as N = number of data points — number of fitted 
parameters —1. Variables that produce the lowest x?ed cor- 
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Figure 7. Rotating cylinder model produced by the best fitting 
parameters. The grey line illustrates a test particle trajectory in 
which vectors show the direction and magnitude of proper mo¬ 
tions and colour indicates the line of sight velocity. 
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Figure 8. Position-velocity diagram for maser features associated 
with the blue lobe of the jet where offsets are measured from 
the jet axis (from Figure [2]l. Also shown is the position-velocity 
profile produced by our rotating cylinder model {solid line). The 
horizontal line at -61 km s” ^ indicates the median velocity of the 
blue lobe masers. 


respond to the best match between model and data. For a 
check of consistency and to investigate the effects of varying 
degrees of freedom on y;?ed considered 3 different condi¬ 
tions when modelling data, which were - to use data from: 

1) all maser featnres. 2) only masers features with measured 
proper motions. 3) only maser features without measured 
proper motions (only position and line of sight velocities). 
The resulting parameters were: 

1) Uout = 45 km s“^, Urot = 22 km s“^, i = 12°: x?ed = 65, 

2) Uout = 45 km s“^, Urot = 19 km s“^, i = 10°: x^ed = 55, 

3) font = 45 km s“^, Urot = 28 km s“^, i = 10°: x?ed = 9. 

All fitting conditions generally agree on a high outflow 
velocity and a shallow inclination angle, with the most de¬ 
viation seen in determination of the rotation velocity. We 
adopt Uout = 45 km s“^, Vrot = 22 km s“^, i = 12° from 
condition 1 results since the fit used all data. 

Values of for parameters using condition 1 are 
shown in Figure [6] where the 3D (wout, Wrot, *) distribution 
of values is represented in ‘slices’ at inclination angles of 
i = 2,12,22°. Modelled and observed maser distributions 
and proper motions are well matched (Figures [2] & [Tjl with 
the redder masers having larger proper motions in a more 
southernly direction and vice versa. Additionally the p-v di¬ 
agram produced by the model shows good agreement with 
the p-v diagram drawn using the maser data (Figure [8]). 

4.4 Overview of jet parameters 

To evaluate the success of the model fitting we calculated the 
scalar sum of 3D motions in the fitting residuals and com¬ 
pared this to the the scalar sum of 3D motions in the data. 
We find that the model can account for 80% of observed 
motions, with 20% remaining as residuals, which we adopt 
as the total model error. The relative fitting success of each 
individual parameter (vout, Vrot and i) can be determined by 
the ‘3D’ shape of the Xred locus in Figure[6l In our fitting it¬ 


erations x?ed increases by 50% its minimum value within 
Auout = ±4, Aurot = ±16 and Ai — ±12, correspond¬ 
ing to expected error contribution ratios of :ef) 

= (1^:4^:3^). The total model error can be separated into 
relative contributors, which allows us to assign error values 
to our best-fit model parameters, which become Uout = 45±2 
km s“^, Wrot = 22 ± 3 km s“^, i = 12 ± 2°. 

Our modelled rotation velocity of Vrat = 22 ± 3 km 
s“^ agrees well wi t h tho se measured bv ICoffev et al.l (l2004h 
and IWoitas et al.l (l2005ll on similar scales. The jet velocity 
gradient is linear and steep, at 1.2 km s“^ mas“^ which 
corresponds to 0.77 km s“^ au“^ at a distance of 1.56 kpc. 
T he very steep gradie nt is consistent with the trend seen 
in IZapata et al.l (l2010l l for velocity gradients to steepen at 
radii closer to the jet. 

The outflow velocity component of the jet is compa¬ 
rable with those frequently reported in the literature and 
the inclination angle is typical fo r a DBSM source, which 
are also known to be jet driven llCaswell fc Phillipj l2008l : 
iMotogi et al.l l2013l . l2015l ~). Regarding the age of the jet, 
if we take a symmetric case where the star stands be¬ 
tween the red and blue lobes the jet length would be r* = 
projected distance/ sin i = 1.5 x 10"^ an, using D = 1.56 kpc. 
Therefore the approximate dynamical age of the maser jet 
is fjet = r*/wout = 1600 years, assuming constant velocity. 
This value will vary between 0 < tjet < (2 x 1600) years 
depending on the position of the star. Note that this is 
not an estimate of the YSO a ge, which is ~ 4000 years 
llDewangan fc Anandar^l2011^ . 

The low uncertainty seen in Uout is echoed by the consis¬ 
tency in the determination of Vout using fitting conditions 1, 
2 and 3. This indicates that LSR velocity information dom¬ 
inates the determination of the outflow velocity. Larger un¬ 
certainty is seen in Vrot this parameter also shows less consis¬ 
tency when comparing the fitting conditions indicating that 
internal motions are important for characterising rotation. 
In the case that a better estimate for the systemic motion of 
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Figure 9. Peak LSR velocity of feature C as a function of time, 
showing acceleration to bluer velocities. 

the source becomes known we may refine the model of water 
masers in S235AB-MIR using more accurate corrections to 
internal motions. 

Although ambiguities may remain in our model due to 
the small number of detected H 2 O masers, it is important 
to note that all models for all combinations of parameters 
and fitting conditions produce results that favour rotation 
in some way, consistently rejecting a ballistic-only interpre¬ 
tation. 

Evidence supporting the presence of a rotating jet is 
itself important as it promotes magnetocentrifugal jet mo¬ 
tions as a viable mechanism for removing angular momen¬ 
tum from an accreting system. We do not quantify this an¬ 
gular momentum using our VLBI data because maser emis¬ 
sion does not scale with mass. Masers have no net charge 
and thus their bulk motions are not directly influenced by 
magnetic fields. We propose that H 2 O masers in S235AB- 
MIR are likely entrained in - and trace the kinematics of 
- a rotating jet of molecular gas. As such, we intend to re¬ 
visit S235AB-MIR using millimeter array observations of jet 
tracers, which should provide suitable quantihcation of the 
angular momentum carried in the jet. 

4.5 Accelerating maser 

A stable velocity increase was observed in feature C, which 
was the only maser feature to exhibit such behaviour. We 
derived a velocity drift rate of 2.92±0.38 km s“^ yr~^, where 
the error value is the standard deviation of a linear fit to the 
data (Figure ID). 

We introduced an additional error term of ±0.63 km s“^ 
corresponding to the half-width of a typical maser spectral 
peak (these are the error bars in Figure[5| to account for the 
possibility that changes in velocity were misinterpreted from 
the changing in dominance of different velocity components 
in a maser feature. Subsequently the measured acceleration 
of feature C is 2.92 ± 1.01 km s“^ yr”^- 

Circular motions give rise to acceleration toward the 
axis of rotation therefore the appearance of an accelerating 
maser is a reasonable finding given the context. The mea¬ 
sured acceleration is consistent with the model prediction. 


within the estimated uncertainty, which supports rotation as 
a possible origin and lends further support to the rotating 
jet interpretation of the velocity gradient in S235AB-MIR. 

The velocity drift rate seen here is similar to that o f 
DBSM maser source G353.273±0.641 (iMotoei et ahllioisl l. 
In their case, slower masers near to the driving source are 
accelerated up to the outflow velocity (120 km s“^), thus 
the masers appear to trace launching motions at the root 
of the jet. In our case, however, the masers seem to be well 
removed from the driving source, supporting an alternative 
expl anation. _ 

iFelli et ^ J 2 OO 6 II also observed a velocity drift in the 
water masers of S235AB using the single-dish Medicina ra¬ 
dio telescope. They find two maser features in the —5 to 
10 km s“^ range which exhibit accelerations of 1 km s“^ 
yr“^. The velocity range creates a near overlap between the 
core veloc i ty an d red lobe masers discussed in this work. 
iFelli et ^ (l2006l l speculate that the velocity drift may have 
been due to the acceleration of masers fro m the central 
YSO up to faster redshift velocities (as seen in iMotogi et al.l 
I 2 OI 5 II . This is certainly feasible given the spatial consistency 
of these maser groups as seen in our VLBI images. Unfor¬ 
tunately these masers were rarely detected in our observa¬ 
tions due to their weak emission. This possible maser outflow 
event would make an excellent target for VLBI observations 
at higher sensitivity. 


5 CONCLUSIONS 

Using VLBI monitoring observations we measured the as- 
trometries of H 2 O masers in S235AB-MIR. Our annual par¬ 
allax of TT — 0.63 ± 0.03 mas corresponds to a distance of 
D = 1.56t°°® kpc, placing the source in the Perseus Arm. 
The source systemic motion was {fiaCosS, jj,s) = (±0.79 ± 
0.12, -2.41 ± 0.14) mas yr"^ 

The S235AB-MIR molecular core drives a maser jet 
however it is unclear if the core contains a single MYSO 
or if it also houses companion lower mass YSOs. As such, 
the identification of the jet progenitor remains pivotal and 
dictates the significance of the VLBI maser data. 

We present evidence that the H 2 O maser burst be¬ 
haviour in S235AB-MIR is caused by the overlapping of 
maser component s alon g the line of sight, as was seen by 
IShimoikura et akl (120051 1 in Orion KL, and propose that sub¬ 
feature scale structures exist which are not fully resolved at 
our current angular resolution. 

Masers trace a young (~ 1600 years) jet which is aligned 
to outflows seen in HCO’^ and C^^S. The maser jet has a 
position angle of PA = 217° and diameter of 15 au, and 
exhibits a well dehned velocity gradient across its axis. The 
jet has a steep velocity gradient of 1.2 km s“^ mas“^ which 
corresponds to 0.77 km s“^ au“^ at a distance of 1.56 kpc. 

3D motions, derived from line-of-sight and sky-plane ve¬ 
locities of 16 water maser features are well described by a 
rotating cylinder model, giving the water maser jet a manner 
comparable to a ‘water spout’. Our best-fit model param¬ 
eters were Uout = 45 ± 2 km s“^, Vrot = 22 ± 3 km s“^, 
i = 12 ± 2°, which are the outflow velocity, tangential rota¬ 
tion velocity and line-of-sight inclination, respectively. Our 
parameters are consistent with other investigations of jet 
rotation at similar scales. These new observations support 
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magnetocentrifugal jet rotation as a possible mechanism of 
dispersing angular momentum from the protostellar disks of 
accreting massive stars. 

One maser feature was found to accelerate steadily at a 
rate of 2.92 ± 1.01 km s“^ yr~^, which was consistent with 
that predicted in our rotating jet model. The velocity drift 
rate is similar to that seen in another DBSM source although 
likely different in origin. 
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APPENDIX A: INVERSE PHASE 
REFERENCING TECHNIQUE FOR VERA 
DATA (MODE: 7 AND 7MM), USING POPS 

As mentioned in Section 2, the i nverse pha s e refe rencing 
(IPR) technique was introduced bv llmai et al.1 (l2012l ~l. Those 
authors provide a discussion of the advantages of IPR and 
the challenges of its implementation in dual-beam VLBI. 
They also provide a link to a webpage hosting working 
scripts that can run IPR (and conventional phase refer¬ 
encing) written in ParselTongue. We introduce here an ap¬ 
proach to IPR VLBI data reduction using AIPS POPS. 

In the IPR approach to reducing VLBI data the phase 
and rate solutions obtained using the maser emission are 
used to calibrate the continuum source data. This practice 
is commonplace for (single beam) VLBI arrays using fast 
switching and a standardised frequency setup, however it is 
more complicated using VERA in VERA7 and VERA7MM 
modes which allocate only a small bandwidth to the maser 
emission and a much wider bandwidth to the continuum 
source. 

VERA dual-beam observations allow simultaneous ob¬ 
servations of maser emission and a continuum source in sepa¬ 
rate telescope beams, henceforth referred to as the ‘A-beam’ 
and ‘B-beam’, respectively. 

With regards to phase calibration the first step is to ap¬ 
ply source specific delay-tracking solutions which are based 
on tropospheric models calibrated usin g GPS measurement s 
of the zenith path delay at each station (IHonma et al.ll2008l) . 
Then we normalise phase gradients in the B-beam baselines 
by solving for the group delay using a bright continuum 
source. These solutions are also applied to the maser data 
since the wide bandwidth of the B-beam provides better 
determination of phase slope than the narrow bandwidth 
of the A-beam. Instrumental phase contributions from the 
dual-beam system were measured using an artificial noise 
source and subsequently removed. 

At this stage all a-priori, station position, instrumental 
and group delay solutions have been applied. What remains 
are mostly atmospheric phase residuals which, at 22 GHz, 
are d ominated by the dynamic troposphere fsee lAsaki et al.l 
l2007f) . 

Phase and rate solutions from the maser emission were 
then obtained with the task FRING in the A-beam at inter¬ 
vals of 1-2 minutes. Solutions were applied to the B-beam 
data - thus phase-referencing the continuum source. Finally 
- the continuum source was imaged, producing a map in 
which the co-ordinates have an origin at the phase tracking 



Figure Al. Flowchart of the inverse phase-referencing procedure 
using VERA7/VERA7MM data. 


center of the B-beam in RA and Dec - with the continuum 
source offset somewhat from the map origin. 

Since the precise co-ordinates of the continuum source 
are well-known its actual offset from the phase tracking cen¬ 
ter should be negligibly small. The measurable offset comes 
from the solutions carried over from the A-beam; the separa¬ 
tion of the maser with respect to the phase tracking center 
in the A-beam. The offset of the continuum source there¬ 
fore gives the absolute co-ordinates of the maser. Fringe- 
rate solutions from the maser were also applied to the A- 
beam, thereby self-calibrating it, which produced high qual¬ 
ity maser maps and revealed many of the weak maser spots 
used in this work. 

The flowchart in Figure Al outlines the reduction 
method in terms of the main AIPS tasks however a step- 
by-step reduction guide, in the format of a POPS log file, 
is hosted at the following webpage: 

http: //milkyway. sci .kagoshima-u. ac . jp/-rossburns88/Scripts .html 
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